Transport in superlattices on single layer graphene. 
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We study transport in undoped graphene in the presence of a superlattice potential both within 
a simple continuum model and using numerical tight-binding calculations. The continuum model 
demonstrates that the conductivity of the system is primarily impacted by the velocity anisotropy 
that the Dirac points of graphene develop due to the potential. For one-dimensional superlattice 
potentials, new Dirac points may be generated, and the resulting conductivities can be approxi- 
mately described by the anisotropic conductivities associated with each Dirac point. Tight-binding 
calculations demonstrate that this simple model is quantitatively correct for a single Dirac point, 
and that it works qualitatively when there are multiple Dirac points. Remarkably, for a two di- 
mensional potential which may be very strong but introduces no anisotropy in the Dirac point, the 
conductivity of the system remains essentially the same as when no external potential is present. 

PACS numbers: 61.46.-w, 73.22.-f, 73.63.-b 
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I. INTRODUCTION 

Graphene is one of the most interesting electronic sys- 
tems to become available in the last few years [l|, 
Graphene is a two-dimensional arrangement of carbon 
atoms in a triangular lattice with two atoms per unit cell. 
In graphene, the electronic low energy properties are gov- 
erned by a massless Dirac Hamiltonian and the carriers 
moving in graphene have very interesting properties: the 
electronic spectrum is linear in the wavevector, and their 
states are chiral with respect the pseudospin defined by 
the two atoms of the crystal unit cell. These properties 
are responsible for exotic effects, such as a half-integer 
quantum Hall effect 0, H| and the Klein paradox - per- 
fect transmission through potential barriers p|. 

The application of electric fields via nano gate geome- 
tries makes it possible to subject the system to poten- 
tials varying on a short length scale. Using these tech- 
niques, recently it has been possible to study experimen- 
tally transport through p-n junctions and p-n-p junc- 
tions in graphene |6h10||. Theoretically, there has also 
been much effort devoted to the study of the spectra and 
the electronic transport through differently doped regions 

EMU 

whose behavior differs from that of conventional 
two-dimensional electron gases. 

A superlattice potential on top of graphene opens the 
possibility of tailoring its band structure and modifying 
its transport properties [l6l - l20j . In particular in the case 
of a one dimensional superlattice potential, the properties 
of the carriers are extremely sensitive to the amplitude 
Vq and period d of the superlattice. For a one dimen- 
sional superlattice, the velocity of the carriers is highly 
anisotropic [2ll423l | and the number of Dirac points at the 
Fermi energy can be altered by varying the product Vod 
[24U261 ] . Moreover, when the potential magnitude of the 
superlattice varies slowly in space, the electronic spectra 
develops a Landau level spectrum [l?} ■ The effect of su- 
perlattice potentials due to external magnetic fields has 



also attracted a great deal of attention [28H31J . 

Several groups have numerically studied electronic 
transport perpendicular to the superlattice barriers [24], 
[H, |32h36| . Starting from the theoretical universal value 



(T = ^^-[LJ], the conductivity increases with the prod- 
uct Vod and develops peaks at the critical values of Vod 
for which new Dirac points emerge [24j|. 

In this work we consider electronic transport in 
graphene in the presence of superlattice potentials that 
are piecewise constant. In the case of one-dimensional 
superlattices we study both transport parallel [Fig. QJa)] 
and perpendicular [Fig. [TJb)] to the barriers. We also 
analyze transport in two-dimensional superlattices [Fig. 
He)]. Analytical expressions for the conductivity are ob- 
tained by describing the carriers with the Dirac Hamil- 
tonian and using the Kubo formula. These are compared 
with numerical results obtained using a tight-binding 
Hamiltonian for graphene in the presence of a superlat- 
tice potential and the Landauer-Biittiker formalism for 
obtaining the electrical conductivity in the presence of 
leads. 

In the case of a one dimensional superlattice, we find 
that, as a function of the product Vod, the conductiv- 
ity parallel to the superlattice barriers, cry, decreases 
quadratically from its value in the absence of the po- 
tential, (To, whereas in the perpendicular direction the 
conductivity <tj_ increases quadratically. The appearance 
of new Dirac points produces peaks in a± and minima in 
cm. For two-dimensional superlattices the conductivity 
depends on the relative values of the product Vod in dif- 
ferent directions. Interestingly, for isotropic superlattice 
potentials, the conductivity is unaffected by the pertur- 
bation and remains at the universal value <to=— tt< Fur- 
ther insight into the character of transport is obtained 
from the channel decomposition of the transmission ma- 
trix. 

This paper is organized as follows. In Section II we 
present the analytical results for the conductivity ob- 
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tained assuming independent anisotropic Dirac points. 
In Section III we present numerical results obtained with 
a microscopic tight-binding Hamiltonian and compare 
with the analytical expressions. Section IV is dedicated 
to the conclusions. 



II. ONE DIMENSIONAL SUPERLATTICE 
POTENTIAL 

1. Preliminaries. 

The electronic structure of an infinitely large flat 
jraphene flake is described by the Dirac Hamiltonian, 



Hq = hvph ■ a 



(1) 



where ?ik is the momentum operator, a are the Pauli 
matrices and v F — 10 6 m/s is the Fermi velocity. The 
two entries of the Dirac Hamiltonian correspond to the 
two carbon atoms in the unit cell in graphene. 

The eigenvalues and eigenvectors of this Hamiltonian 

are ek, s =sv F hk and \s, k >= £-7=- ^ se i9(k) 1 > where s — 

— 1 and s = 1 describe the occupied and empty bands 
respectively. In the previous expressions 0(k) is the angle 
of the vector k with respect to the k x direction. 

2. Superlattice band structure. 



where V- 



Vpd 



- 2Hv . The group velocity of the state is 
anisotropically renormalized, and has a strong depen- 
dence on the direction of the wave vector k[2l|. At 
the Dirac point and for directions along the superlat- 
tice axis the velocity of the carriers is unaffected by the 
potentials, v x =v F . However the group velocity along the 
direction perpendicular to the superlattice direction is 
strongly renormalized and takes the form 



v F - 



\sm(V)\ 



V 



(3) 



Whenever the superlattice parameters satisfy the con- 
dition, 



Vpd 
hvp 



2nj j = 1,2,3, ... 



(4) 



the group velocity in the y direction vanishes and a 
new pair of Dirac points emerges from the original 
Dirac point, moving in opposite direction along the k y - 
direction[2J, [25[. Near the new Dirac points and at low 
energy the dispersion is also linear and anisotropic. For 
the j-th pair of new Dirac points the velocity in the x 
and y directions have the expressions (26j], 
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(5) 
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FIG. 1: (Color online) Schematic representation of the su- 
perlattices used showing the axis selection. The system is 
infinite along the ^-direction and has a finite length along 
the {/-direction. The dashed patterns on each side of the 
y-direction indicate the leads for the Landauer conductance 
calculations. The superlattice barriers can be parallel (a) or 
perpendicular (b) to the direction of transport. We also con- 
sider a chessboard-like two dimensional superlattice potential 



in which cL 



d x (c). 



We consider a one-dimensional Kronig-Penney super- 
lattice along the x-direction (see Fig. Ufa)). The 
period of the potential is d, V(x) — V(x + d) and 
V (x)=Vq sgn(x) for \x\ < d/2. For this potential it is 
possible to find an analytical expression for the band 
structure [ljl [26j], that in the limit of small wave vec- 
tor and energies takes the form 



e(k) = ±Hv F k\ 



sin 2 (f ) 
V 2 



1/2 



(2) 



3. Electrical conductivity. 

The conductivity in the collisionless limit has the ex- 
pression (371. |38|| 



e_ ^ h,s'~fk,s I < s, k\v^\s\ k > - 

k,s,s 



-£k,s £k„ 



-£k,s 



-huj — i8 



where s' and s are band indices, /k, s is the Fermi distri- 
bution function for the states \s,k >, is the velocity 
operator in the /t direction and 5 is a positive infinitesi- 
mal constant. The conductivity contains a factor 5 s .g u =4, 
which takes into account the spin and valley degeneracy. 
In the case of a single Dirac point with anisotropic ve- 
locities v x and v y , expressed with a Dirac Hamiltonian of 
the form 

HA = h(V X k x a X + VykyCTy), 

one may show that the conductivity parallel and perpen- 
dicular to the potential barriers of the superlattice may 
be written in the form 



af[(w = 0) = 
a° ± (uj = 0) = 



"(TO = (TO" 



\sm(V)\ 



„0 



(Tq = (Tq- 



V 
V 



sm(V)\ 



(7) 



3 



with do the conductivity of an isotropic Dirac Hamilto- 
nian. The value of o~q depends on the order in which the 
zero frequency, zero temperature and vanishing "smear- 
ing parameter" 5 [38| limits are taken(38l l39j. However 
the form of the velocity rescaling of the conductivity is 
independent of the order in which the limits are taken. 

In the case of several Dirac points in the spectrum, we 
assume that each of the points contributes to the con- 
ductivity in parallel and using Eq. ([5]), the conductivity 
takes the form, 



C|| = Co 



<7j_ = cr 



V 
V 



+ *E 



V 2 (njf 



2 E- 



|sin(U)| ^_V*-(irj) 



(8) 



where j max =hxteger(—) indicates the number of Dirac 
point pairs induced by the superlattice. From this ex- 
pression we see that for small potentials the conductiv- 
ity perpendicular to the superlattice barriers increases 
quadratically with V^d, and each time a new pair of Dirac 
points emerges the conductivity exhibits a peak. In the 
direction parallel to the barriers, the conductivity de- 
creases quadratically with Vod and dips when new Dirac 
points emerge. 

We remark that in obtaining Eq. ©, we have as- 
sumed that each Dirac point contributes as an indepen- 
dent channel to the conductivity and that near each Dirac 
point the dispersion relation is linear over a wide range 
of the reciprocal space. 

4- Mode dependent transmission. 

The conductivity of a system governed by the Dirac 
equation with anisotropic velocities, H = h(v x k x a x + 
VykyGy), can be also obtained by calculating the trans- 
mission probability of modes confined in a stripe of width 
W and length L connected to heavily doped contacts [4p|— 
|42| . For transport along the i-direction, the transmission 
probability for a transverse mode has the form 



T n (x) 



1 



cosh 2 (^q n L) 



(9) 



where the transverse momentum q n depends on the 
details of the precise boundary condition of the strip 
(iol I43I ] . For wide enough strips the conductivity of the 
system is independent of the boundary conditions and is 
found by summing over the modes, 



Le 2 ^ e 2 2L 
9s9v— 7 — T n\X, 



W h 



dq 



4 e 2 v x 
7T h v„ 



for W » L . 



h 7T cosh 2 (^qL) 

(10) 



The conductivity in the y direction is obtained by inter- 
changing x and y in the last equation. The condition for 



the existence of a well defined -size independent- conduc- 
tivity is the dependence of the transmission probability 
on the product qL (Eq. ((9|)) and the linear dispersion of 
the carriers. The condition W ^> L allows the sum the 
transmissions over the modes to be written as an integral 
over q in Eq. (fl0|l . 



5. Two dimensional superlattice potential. 

We consider a two dimensional superlattice potential 
on top a graphene sheet [as in Fig. QJc)]. In second order 
perturbation theory the group velocity of quasiparticlcs 
with momentum k has the form 12 ill 



2|£/(G)| 2 . a 



'k,G 



(11) 



where G and U(G) are the reciprocal lattice vectors and 
the corresponding Fourier component of the external po- 
tential and 6>k.G is the angle between G and k. Using 
the same approximation as in the previous subsection 
the conductivity in the x-direction takes the form 



cr - 



h 2 v 2 F -E G ^\U(G)\ 2 j^ 

W-Eg^I^gOPtS* 



(12) 



The conductivity in the y-direction is obtained by in- 
terchanging G x and G y in this expression. The strik- 
ing result of Eq. fj 12|) is that for symmetric superlat- 
tice potentials the conductivity in the x and y direc- 
tions are equal and take the value of pristine graphene, 
o~ xx =o- yy =o-Q=^^ . The expression Eq. (fTT|) has been 
obtained in second order perturbation theory and it is 
a good approximation provided that the superlattice po- 
tential does not induce new Dirac points. We expect that 
Eq. ([T2"T) will be valid in the same regime. 



III. NUMERICAL CALCULATIONS. 

In order to compute numerically the transport prop- 
erties we describe the electronic states of a defect free 
graphene layer using the tight-binding approximation, 



H 



-t, 



E 

<ij> 



•E^ici 



(13) 



where t g = 2hvF/3ciQ denotes the hopping element be- 
tween nearest carbon atoms on the hexagonal lattice, do 
is the smallest carbon-carbon distance and Vi is the po- 
tential applied to the lattice. The spin degree of freedom 
has been omitted due to degeneracy. 

In order to analyze the different transport situations 
depicted in Fig. Q]we assume that the central region is a 
nanoribbon with armchair edges along the x-direction as 
depicted in Fig. [21 The nanoribbon is constructed by re- 
peating a unit cell composed of four atoms N times along 
the y-direction and M times along the ^-direction. Thus, 
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FIG. 3: (Color online) In the left panels we plot, as function of 
qd, the transmission T(q) per spin channel for a superlattice of 
period d — 42ao and amplitudes V = or" = an d ^ ~ 0.7 
(top left panel) and for a period d = 54ao and amplitudes 
V — 4.3 and V = 5.4 (bottom left panel). In the right panels 
we plot the distribution of the eigenvalues of the transmission 
matrix for the different values of V. The length of the stripe 
is L = 100V3o . 



FIG. 2: (Color online) Schematic representation of the super- 
lattice used in the tight-binding calculations showing the axis 
selection. The system is infinite along the x-direction and 
has a finite length L = y/3Nao along the {/-direction. The 
superlattice has a vertical period d = 3Moo in which periodic 
boundary conditions are imposed. 



the length of the graphene stripe is L = N\/3aQ. For 
describing the W 3> L limit we impose periodic bound- 
ary conditions in the transversal direction x and define 
q € [— n/dsc, ^/dsc] as the corresponding wave vector, 
with dsc = 3Mao being the vertical length of the super- 
cell. 

We connect the armchair edges of the nanoribbon 
to heavily doped graphene leads thus maintaining the 
graphene sublattice structure at the edges |44Tj48| . The 
corresponding self-energies on the graphene sites at the 
layer edges are approximated by a AM x AM matrix with 



elements T 



L.R 



ij.afi — ^ijla'p 1 ' wuere a, j8 = 1, . . . , 4 label the 
atomic sites within the unit cell and i,j — 1, . . . , M label 
the unit cells in the superlattice. Following the geometry 
depicted in Fig. [2 the elements of the self-energy matrix 
are explicitly defined as 7^ 2 = 7^3 ==Ju = 744 = «v3/2 
and 723 = 732 = 7/4 = 7I1 = — l/2[47j. Thus, we calcu- 
late the transmission at zero energy, T(q), as 



T(q) = 4Tr f L G r LR (E = 0, q)f R G a RL (E = 0, q) , (14) 



mission matrix iH, where i = 2yT lG t lr (E, q) y T R . 
From these eigenvalues one can determine the probabil- 
ity distribution P(t) = J2 a q ~ T a{l)) an d the Fano 



factor 



F 



iM 

£ Er a (q)(l-r a (q)) 

a— 1 q 

' 4M 

EEr«(9) 

a—1 q 



(15) 



By integrating the transmission we compute the conduc- 
tance of the system G = (Ae 2 /h) J dqTr , where both 
the spin and valley degeneracies have been taken into 
account. The resulting conductivity, within the limit 
W S> L, is obtained by multiplying by the geometrical 
factor L. 



A. Transport parallel to the superlattice barriers. 

For studying the transport parallel to the superlat- 
tice, we consider a periodic one-dimensional potential 
along the ^-direction within the previous geometry as is 
schematically depicted in Fig. QJa) . The one-dimensional 
superlattice potential, Vi, has the piecewise constant 
form, 



where G r £ R RL (E, q) are the AM x AM retarded and ad- 
vanced Green functions between the edges of the layer. 
Furthermore, for analyzing the transmission distribution 
it is useful to determine the eigenvalues r Q (q) of the trans- 



Vi = 



V 
-Vo 



< a* < § 

1 < Xi < d 



(16) 



where d 



>-sc 



3Moq is the period of the potential. 
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In Fig. [3] we plot the transmission T(q) as function 
of the product qd for a superlattice of period d = 42ao 
and amplitudes Vod = and Vod — IAHvf in the top 
left panel and for a period d — 54ao and amplitudes 
Vod = 8.6hvp and Vod — !0.8hvF in the bottom left 
panel. The horizontal length of the graphene layer is 
L = 100\/3ao- We also plot in the right panels of Fig. 
12] the distribution of the eigenvalues of the transmission 
matrix. 

In Fig. [4] we plot, as function of Vyd, the conductiv- 
ity and the Fano factor obtained for a system of length 
L=100V3ao a nd for different values of the superlattice 
period, d. 

We first discuss the case of potential barriers in the 
range V < V c = [top panels of Fig. [3J. For these 

superlattices the original Dirac points are the only active 
transport channels. As a function of q the transmission 
is peaked at q=0, and the width of the peak diminishes 
when Vod increases. The transmission fits very well to the 
functional form [see Eq. ©] T(q) = 2/ cosh 2 ( — qL), 

where the factor 2 accounts for the valley degeneracy 
and V = ^f. The corresponding distribution of the 
eigenvalues of the transmission matrix has the form 
P(t) ~ l/Vyl — t indicating the pseudo-diffusive char- 
acter of the transport in this range of potentials. The 
conductivity obtained by integrating the transmission 
is well-defined and, in this range of Vod, has the form 



(TO 



|sin(y)| 
V 



[see Fig. [4]. The Fano factor in this 
range of potentials is 1 /3 in agreement with the pseudo- 
diffusive character of transport. We thus conclude that 
in the range of parameters Vod < 2'KhvF the transport 
is pseudo-diffusive, the conductivity only depends on the 

product Vod and has the form a® = o~o ^ blD ~ v ^ . 

For normalized barrier heights V>d larger than 2nhvF 
two new Dirac points per valley appear [241 [25| . These 
new Dirac points are new transmission channels in the 
system, that for transport parallel to the superlattice bar- 
riers are superimposed in reciprocal space upon the orig- 
inal Dirac points. The resultant transmission exhibits 
a wider distribution in reciprocal space [see bottom left 
panel of Fig. [3J. The width of the transmission can 
reach the edges of the reduced Brillouin zone ±ir/d for 
small values of L/d. The corresponding distribution of 
the eigenvalues of the transmission matrix is a superposi- 
tion of the distribution of each mode, and the correspond- 
ing Fano factor is different than 1/3. The conductivity 
should be independent on the system size. We find that 
the value of L where the conductivity is well defined de- 
pends on d and coincides with the value of L in which 
the transmission is non zero at the edges of the reduced 
Brillouin zone. In Fig. @]we see that the general trend of 
the conductivity for values of Vod larger than 2'KhvF is 
qualitatively described by the continuum model, Eq. . 
However the analytical model neglects some effects such 
as the coupling between the modes or the deviation from 
linear dispersion, so that in this range of superlattice pa- 
rameters the conductivity depends separately on Vq and 
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FIG. 4: (Color online) Transport parallel to the superlattice 
barriers. Top (bottom) panel shows the conductivity (Fano 
factor) for a graphene sheet with L = 500\/3ao and superlat- 
tice period d — 24ao (solid blue line) and L — 200\/3ao and 
superlattice period d = 36ao (dashed red line) as a function of 
the normalized barrier height Vod. Dotted line corresponds to 
the conductivity obtained in the continuum model assuming 
independent transport channels, Eq. in the top panel and 
to the pseudo-diffusive value F = 1/3 in the bottom panel. 



d. The coupling between the modes also leads to a Fano 
factor with a value larger than 1/3, and the transport is 
not pseudo-diffusive. 



Transport perpendicular to the superlattice 
barriers. 



In this section we consider a potential in the y-direction 
and study the transport in the same direction, i.e. per- 
pendicular to the superlattice barriers [see Fig. QJb)]. 
Following the same geometry as in the previous section 
(see details in Fig. [2]), we define a one-dimensional piece- 
wise potential along the y-direction as 



Vi 



Vo 
-Vo 



0<yt 

i<yi 



< 



d 

2 

< d 



(17) 



where d = 2n\/3ao is the period of the potential. 

In the left column of Fig. [5] we plot the transmission 
T(q) as a function of qd$c f° r a superlattice with period 
d = 38.1ao and amplitudes Vod = 0, Vod — 5hvp (top 
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FIG. 5: (Color online) In the left panels we plot, as func- 
tion of qdsc, the transmission T(q) per spin channel for a 
superlattice of period d — 38.1ao and normalized amplitudes 
V = ^ = 0,V = 2.5 (top left panel), V = 3.3 and V = 3.6 
(bottom left panel). In the right panels we plot the distri- 
bution of the eigenvalues of the transmission matrix for the 
different values of V. The length of the stripe is L = 500\/3ao. 



left panel), Vod — 6.6hvp and Vod — 7.2hvp (bottom left 
panel). The horizontal length of the graphene strip is 
L = SOOV^ao- In the right column of Fig. [5] we plot 
the corresponding distribution of the eigenvalues of the 
transmission matrix. 

In the top panel of Fig. [5] we show, as function of Vod, 
the conductivity for horizontal periods of d = 34.6ao and 
d = 76.2ao, for a graphene sheet of length L=500\/3ao- 
In the bottom panel of Fig. [6] we plot the Fano factor for 
the same two values of the period of the superlattice. 

In the range of potential barriers before the creation 
of new Dirac points, i.e. Vq < V c , the behavior of the 
transmission is exactly the inverse of the previous case. 
The contribution to the transmission from each valley 
is superimposed as a sharp peak at q = 0. However, 
contrary to the previous result, the width of the peak 
increases with the product Vod. Following Eq. ©, the 

transmission is fitted to T(q) = 2/ cosh 2 ( ■ sn l^ v \ gX). Sub- 
sequently, the distribution of the eigenvalues is that of 
pseudo-diffusive transport. On the other hand, when 
Vo > 27T hvp, a pair of Dirac points is created for each val- 
ley. In the bottom panel of Fig. [S]we show how these new 
peaks split from the original ones until there are three al- 
most independent contributions to the transmission. In 
this later case the distribution of eigenvalues for each 
mode returns to a form of the type P(t) ~ 1/r-^l — r, 
indicative of pseudo-diffusive behavior. Before the new 
Dirac points are completely separated from the original 
ones, the coupling between modes produces a deviation 
from the pseudo-diffusive transport. 

The behavior of the conductivity perpendicular to the 
barriers is completely different than the parallel case. 
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FIG. 6: (Color online) Transport perpendicular to the super- 
lattice barriers. Top (bottom) panel shows the conductivity 
(Fano factor) for a graphene sheet with L — 500\/3ao and 
superlattice periods d — 34.6ao (solid blue line) d — 74.2ao 
(dashed red line) as a function of the normalized barrier height 
Vod. Dotted line corresponds to the conductivity obtained in 
the continuous model assuming independent transport chan- 
nels, Eq. iJHJ, in the top panel and to the pseudo-diffusive 
value F = 1/3 in the bottom panel. 



The perpendicular conductivity presents peaks at the val- 
ues of the normalized potential height where new Dirac 
points appear. The numerical calculated conductivity 
agrees very well with the analytical one, Eq. (j8]), even 
for values of Vod > 2irhvp. The Fano factor has the value 
1/3 for all values of Vod except near the values of Vod for 
which new Dirac appears. This indicates that, in this 
geometry, the Dirac points are weakly coupled and the 
approach of Section II for the conductivity is appropriate. 



C. Transport in a two dimensional superlattice 

One of the more striking results presented in Section II 
is that the conductivity of graphene in the presence of a 
symmetric two dimensional superlattice potential is — j- 
independent of the period and the height of the potential 
barriers. In order to check this result we have built a 
chessboard-like potential combining piecewise potentials 
in the x, Eq. (|16[) . and y, Eq. (|17[) . directions in a way 
in which a potential barrier is always followed by a well 
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FIG. 7: (Color online) In the top (bottom) panel we plot 
the conductivity (Fano factor) as a function of the normal- 
ized barrier height Vod y for a graphene sheet in presence of 
a two dimensional superlattice with a fixed vertical period 
of d x = 48ao and different horizontal periods d y = 45ao 
(blue solid line) and d y — 48.5ao (red dashed line). The dot- 
ted line corresponds to the conductivity of pristine graphene 
(<to) in the top panel and to the pseudo-diffusive limit (1/3) 
in the bottom panel. The length of the graphene sheet is 
L — 200\/3fflo. Inset: the conductivity as a function of the 
normalized barrier height in the proximity of the critical po- 
tential V c in which a new pair of Dirac points is created. 



along each direction (see Fig. Die)). The length of the 
period in the x and y directions is d x and d y respectively. 
Because the underlying triangular lattice of graphene, the 
period in both directions cannot be exactly equal. 

The top panel of Fig. [7] shows the conductivity as 
a function of the potential height Vq for a graphene 
layer with L = 200v / 3ao and a fixed vertical period of 
d x = 48ao- We plot the conductivity for two differ- 
ent horizontal periods d y = 45ao and d y — 48.5ao. We 
compare these results with the isotropic conductivity of 
graphene <r = 

A remarkable result is that the conductivity in this 
potential remains almost constant in the range Vb > V c 
where a new pair of Dirac points is created in the previ- 
ously studied cases. Thus, in this range of potential bar- 
riers, Eq. (1121) obtained in second order perturbation the- 
ory remains a good approximation according to the tight- 
binding results. Furthermore, the pseudo-diffusive be- 



havior of transport is maintained for a large range of the 
potential barriers. In the bottom panel of Fig. [7] we show 
how the Fano factor is stable around the pseudo-diffusive 
value of 1/3 while Vod y < AttHvf- When Vod y ~ AttHvf, 
which for the previous potentials corresponded to the cre- 
ation of the second pair of Dirac points, the conductivity 
deviates from <7q, the Fano factor increases and trans- 
port is no longer pseudo-diffusive. The approximation of 
weakly coupled Dirac points is then no longer applicable. 

The small deviations from the conductivity of pristine 
graphene that occurs when Vod y ~ 2ttHvf can be more 
clearly appreciated in the inset of Fig. [7] Due to the 
geometry of the graphene layer, the period in both direc- 
tions is never exactly the same. This affects the validity 
of Eq. (|12p to a small degree. When d y < d x the con- 
ductivity slightly increases from <7o, presenting a positive 
slope, while if d y > d x the effect is the opposite. When 
the difference between both periods becomes larger the 
conductivity continuously evolves into the corresponding 
case of the previous sections (Figs. I4l6p . 



IV. CONCLUSIONS 



Superlattice potentials generically induce anisotropy 
in the dispersions near the Dirac points in graphene, 
and under certain circumstances may induce extra Dirac 
points at zero energy. In this work we demonstrated that 
when the Fermi energy passes through a spectrum with 
a single anisotropic Dirac point, the resulting conductiv- 
ity can be expressed in a very simple way in terms of 
the velocities along the two principle directions of the 
anisotropy, and the conductivity for the corresponding 
isotropic Dirac point. The result can be generalized to 
the case of several Dirac points when they are sufficiently 
separated in momentum space so that a conductivity ex- 
pressed as a sum over those of independent Dirac points 
is sensible. For a two-dimensional superlattice which in- 
duces little anisotropy in the spectrum, a remarkable re- 
sult is that the conductivity is essentially unchanged from 
the result for pristine graphene, even if the velocity renor- 
malization is quite large. 

Numerical tight-binding calculations generally confirm 
this simple picture. In particular one finds the conductiv- 
ity parallel and perpendicular to the superlattice barriers 
for a one-dimensional potential evolve in opposite direc- 
tions with increasing Vod, and that for a spectrum in 
which no new Dirac points have been generated there is 
quantitative agreement with the simple analytical model. 
As new Dirac points are introduced into the spectrum one 
finds dips in cry and peaks in a± as expected, although the 
results are less quantitatively described by the continuum 
model, presumably because the wavefunctions cannot be 
uniquely associated with single Dirac points. Deviations 
of the Fano factor from pseudo-diffusive behavior confirm 
this interpretation. 
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These studies suggest that more complicated poten- 
tials could also yield behaviors in the conductance with 
simple interpretations. For example, a modulated super- 
lattice potential yields a Landau level spectrum [27|, for 
which (Tii may have behavior reminiscent of edge state 
transport 49]. It is also interesting to speculate that 
for isotropic superlattice potentials, one may sufficiently 
slow the electron velocity so that electron-electron inter- 
action effects become important [5(1 [5l[ . We leave these 
questions for future research. 
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